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Abstract 

An astonishingly simple analytical frequency approximation formula 
for a class of nonlinear oscillators with large amplitudes is derived and 
applied to various example systems yielding useful quick estimates. 


1 Introduction 

In addition to established methods like Harmonic Balance, Krylov Bogoliubov 
or Lindsted Poincare m many new approaches for approximating the limit 
cycle frequencies of strongly nonlinear oscillators have been introduced in recent 
years, e.g. the Energy Balance method m, the Hamiltonian Approach m, the 
Variational Iteration method m , the Amplitude frequency formulation m or 
the Newton Harmonic Balance Method m and other methods P5l liq . These 
new methods have been successfully applied to various systems (see e. g. |U 
i m h nu 0 on). Here, we present an extremely simple straightforwardly 
applicable analytical frequency approximation formula which yields satisfactory 
results for a variety of systems and parameter ranges with a minimum of effort. 
It is particularly useful in cases where the potential (or antiderivative of the 
force function) cannot be written in an analytically closed form. 

As a first example, let us consider the Duffing oscillator 


x + ax + ex 3 = 0 (1) 

with the initial conditions x(t = 0) = A, x(t = 0) = 0. The ansatz 

x(t) = Acos(ujt) (2) 

satisfies the initial conditions and becomes an exact solution in the linear case 
e = 0 if cu 2 = a. Inserting our ansatz (J2J into the differential equation JTJ we 
obtain 

— v 2 Acos(ujt) + oAcos(wf) + eA 3 cos 3 {cot) = 0 (3) 
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where the cubic cosine function can alternatively be written as cos 3 (wt) = 
\ (3 cos(wf) +cos(3wf)) pQ. Due to the term proportional to cos(3 ut) our ansatz 
x = Acos(<jot) cannot be an exact solution of equation (JT]) . We seek an approx¬ 
imate solution for the frequency w by means of a colocation method, i.e. by 
evaluating equation (0 at some time t £ [0,T/4] where T = 2it/u>, similar to 
the procedure used in m in the context of the Energy Balance method. In m 
He used an analogous approach in combination with a Galerkin method rather 
than colocation. 

We want to choose our colocation time t such that the influence of the 
cos(3wt)-term is small. We therefore evaluate the differential equation where 
cos(3 uit) = 0 for the first time which leads to the condition 3wf = 7r/2 or 



(4) 


The colocation point wt = 7r/6 was also successfully used, on a purely penomeno- 
logical basis, in the context of He’s Amplitude Frequency formulation mm (see 
also next section and the appendix). Inserting condition 0 into equation 0 
we obtain _ 

to = a + eA 2 cos 2 ^ a + A 2 (5) 

with cos(7r/6) = x/3/2. This approximate result coincides with other ap¬ 
proaches like first order Harmonic Balance m , first orders of He’s Energy 
Balance method M and his Hamiltonian approach m as well as other meth¬ 
ods [m na ini [23] . 

The energy of the Duffing oscillator 0 is given by E = ^x 2 + V(x) with the 
potential V(x ) = jx 2 + |x 4 . This leads to the expression w exact = 27 t/T exact 

where [ 22 ] T exact = f* A sj dx = N K ( r "’- m = 2 (*+eA*) and K ( m ) 
is the complete elliptic integral of the first kind. 

Comparing the approximation (0 with the exact frequency w e xact one can 
show m that the relative error of u> is always less than 7.6% even in the extreme 
large amplitude limit eA 2 —> 00. 

The rest of the artricle is organized as folows: In section 0 generalizing 
the calculation from the introductory example a simple frequency approxima¬ 
tion formula for nonlinear oscillators with antisymmetric force terms f(x) which 
have a Taylor expansion at x(0) = A is derived. In section [3] the application 
of the formula is illustrated by means of several examples, including nonlinear 
oscillators described by discontinuous force functions and force functions with¬ 
out closed form potential energies. In section [4] the method is generalized to 
asymmetric force functions. The main results are summarized in section 0 In 
appendix [0 we demonstrate that the frequency approximation formula derived 
in this paper is mathematically equivalent to a certain phenomenological variety 
of He’s amplitude frequency formulation. 
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2 Simple frequency approximation formula 

Now we consider a more general nonlinear oscillator of the type 


x + f(x) = 0 ( 6 ) 

with the initial conditions x(t = 0) = A, x{t = 0) = 0 and where f(x) is 
antisymmetric in x , i.e. f{—x) = —f(x) and differentiable in x ^ 0. Due to the 
symmetry properties of the system the Fourier expansion 

OO 

*«) = £ b 2 k-i cos((2fc - l)ut) (7) 

fc=l 

contains only odd multiples of wt 111: (see also the discussion in [27j where the 
same class of systems is considered). Thus the leading and next to leading order 
terms are cos(wt) and cos(3 cot) respectively. 

Let us further assume that the force function f{x) can be Taylor-expanded 
in the vicinity of the initial amplitude x(0) = A, f(x) = f(A ) + f'(A)(x — 
A) + O ((x — A) 2 ). The Taylor expansion of x(t) can be written as x(t) = 
A — \a 2 t 2 + jidit 4 + 0(t 6 ). Inserting these expansions in ([ 6 ]) yields a 2 = f(A), 
04 = f'(A) leading to x(t) = A — | f(A)t 2 + jff'(A)f(A)t 4 + 0(t 6 ). Comparing 
this expression with the Taylor expansion of the Fourier representation in 0 
we obtain A = J2kLi b 2 k-i, f(A) = uj 2 ~ l) 2 ^ 2 fc-i and f(A)f(A) = 

w 4 E^Li( 2 fc- l) 4 ^ 2 fc-i- Convergence of the second (resp. third) of these series 
implies limfc_>, 00 ( 2 /c — l) 2 & 2 fc-i = 0 (resp. linp c _ ) . 00 ( 2 fc — l) 4 & 2 fc-i = 0 ) such that 
b 2 k-i 0 faster than l/( 2 fc — l ) 2 (resp. b 2 k-i —> 0 faster than 1 /( 2 k — l) 4 ) in 
the limit k —> 00 . This motivates approximations based on a truncation of the 
Fourier series 0 after a finite number of terms. 

As in the introductory example we insert the ansatz x(t) = Acos(cot) into 
our differential equation 0 arriving at 

— co 2 Acos(ujt) + f(Acos(ut)) = 0 . (8) 

In analogy to the introductory example we colocate at cot = 7 r/ 6 , where the 
next to leading order terms proportional to cos(3 uit) are zero, which leads to 
the simple approximation formula 


U! = 


N 



(9) 


In appendix El we show that this approximation formula is mathematically 
equivalent to a certain variety of He’s amplitude frequency formulation based 
on phenomenologically chosen parameters. 
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A CJRK [6 ^approx Error (%) 


1 

2.2798 

2.3049 

1.1010 

5 

2.7318 

2.7500 

0.6662 

10 

3.2057 

3.2210 

0.4773 

100 

7.7762 

7.8102 

0.4372 

1000 

23.7999 

23.8170 

0.0718 


Table 1: The approximate frequency w appr ox from (fllT) for the cubic-quintic 
oscillator is compared with frequencies wrk from numerically exact Runge Kutta 
calculations [6j for different values of A and a = 1, e = 5, A = 1. 

3 Example applications 

3.1 Example 1 

The cubic quintic oscillator with the force function 

f(x) = ax + ex 3 + Acr 5 (10) 

reduces to the Duffing Oscillator © in the limit case A = 0. The simple ap¬ 
proximation formula © yields the frequency 

uj = ^a + e^A 2 + A^A 4 . (11) 

A comparison with numerically exact frequencies for different values of A (table 
UJ reveals a good agreement. 

3.2 Example 2 

The fractional strongly nonlinear oscillator described by 

/Or) = x 1 ' 3 , (12) 

has been considered in several articles Gramms]- From equation © we obtain 
the approximate frequency 

1/6 

w=(-j A~ 1/3 « 1.0491A“ 1/3 (13) 

which coincides with the first order Harmonic Balance result IB]. A comparison 
with the exact frequency uj ex = 1.070451A -1 / 3 G3 reveals an error of approxi¬ 
mately 2.0%. 

3.3 Example 3 

Next we consider the strongly nonlinear oscillator with 

f{x)=x ~ 1 (14) 
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analyzed in mmm- Formula © yields the approximation 


w = -^A- 1 « t.1547^- 1 (15) 

V3 

coinciding again with the first order Harmonic Balance result [20] and results 
from He’s homotopy perturbation method [T2j . The exact frequency reads 
w ex = ss 1.2533141A -1 [20] . The resulting error of 7.9% is accept¬ 

able considering the simplicity of our approach. 

3.4 Example 4 

An oscillator with the force function 

fix) = x - A X (16) 

V1 + x z 

was used in HEHElj to model the dynamics of a mass attached to a stretched 
wire. Using formula © we obtain the approximate frequency 

w = J 1- ; A (17) 

V v/1 + (3/4)A 2 

coinciding with the result obtained in [5] using a Harmonic Balance approach, 
which was shown to be in good agreement with the exact frequencies yielding 
relative errors no larger than 2.7% for 0 < A < 1. 

3.5 Example 5 

For the nonlinear oscillator with the force function 

/(x) = sin(x 3 ), 0 < x 3 < pi/2 (18) 

the potential energy cannot be written in an analytically closed form. The 
approximation © reads 


w = 


sin(3 3 / 2 A 3 /8) 




(19) 


A comparison with Runge Kutta calculations (tabled shows the validity of the 
approximation m for values 0 < x < 1. 


3.6 Example 6 

Next we consider the force function 

/(x) = exp(/3-y/[x|)sign(x) (20) 
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A UJRK c^approx Error (%) 


0.1 

0.0847 

0.0866 

2.2 

0.5 

0.423 

0.433 

2.4 

0.75 

0.630 

0.645 

2.5 

1 

0.806 

0.836 

3.7 


Table 2: The approximate frequency u) appr0 x from (fTTH) for Example 5 is com¬ 
pared with results of a Runge Kutta integration ujrk for different values of 
A. 


P = — 1 : P = 1: 


A 

urk 

^approx 

Error (%) 

A 

IaJRK 

^approx 

Error (%) 

0.1 

3.0448 

2.9331 

3.7 

0.1 

4.0525 

3.9367 

3.0 

1 

0.7073 

0.6748 

4.6 

1 

1.7470 

1.7112 

2.1 

5 

0.1828 

0.1698 

6.3 

5 

1.3713 

1.3602 

0.81 

10 

0.0845 

0.0780 

7.6 

10 

1.4813 

1.4800 

0.088 

Table 3: 

The approximate frequency w a ppr0 x from (HID) for 

Example 6 

is com- 


pared with results of a Runge Kutta integration ujrk for different values of A 
with P = — 1 and p = 1. 


which is discontinuous at x = 0. Furthermore it potential energy cvannot be 
written in terms of elementary functions if p ^ 0 (where sign(a;) = 1 if x > 0, 
sign(a;) = — 1 if x < 0). Formula © leads to the approximation 


w = 



N 


4a 


( 21 ) 


For P = 0. comparison with the exact solution m Vexact = yields an 

error of approximately 3.3 %. Comparison with Runge Kutta calculations for 
P = ±1 (see table O also reveals a decent agreement. 


4 Generalization to asymmetric force functions 

Now we drop the requirement of antisymmentry and consider a more general os¬ 
cillator with rest position x = 0, /( x) continuous if x ^ 0 and Taylor-expandable 
about x = ±A. According to formula ([9|) we define the frequencies 


w± = 



( 22 ) 


where x± denote the right and left turning points. We consider the initial 
conditions #(0) = £+, i(0) = 0. Following Hu [T5] we decompose the half 
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A 

Texact [IS] 

Tapprox 

Error (%) 

0.1 

6.3116 

6.3122 

9.506 e-3 

0.2 

6.4414 

6.4140 

4.254 e-1 

0.3 

6.6294 

6.6354 

9.051 e-2 

0.4 

7.1246 

7.1318 

1.011 e-1 

0.45 

7.7065 

7.6967 

1.272 e-1 

0.49 

9.2080 

8.9909 

2.358 e+0 

Table 4: 

. The Period Tapprox 

from (1231) with right turning point x+ 


turning point (1251) for the asymmetric oscillator with f(x) = x + x 2 is compared 
with exact results T exa ct m for different values of A. 


period T/2 into two parts: the elapsed time between x = x+ and x = 0 is 
approximately given by > the time between x = 0 and x=£_ by ■ Thus 



is an approximate value for the period and 



(23) 


(24) 


for the frequency of the oscillator. 

As an example we consider the asymmetric oscillator with f(x) = x + x 2 
with initial conditions x(0) = x + = A > 0, i(0) = 0. Energy conservation 
E = ^ x ± + ^ x ± leads to [T5j 


x _ 




(25) 


For the existence of a periodic solution the discriminant in (1251) must be pos¬ 
itive which leads to the condition 0 < A < 1/2. In table [4] we compare the 
approximation (1231) with exact results m for several values of A, obtaining a 
good agreement. 


5 Conclusion 

A simple frequency approximation formula for a class of strongly nonlinear 
oscillators with antisymmetric position-dependent force terms was derived and 
applied to several example systems. This derivation provides a justification for 
a hitherto only phenomenologically applied variety of the amplitude frequency 
formulation which is shown to be mathematically equivalent to the formula 
presented here. The formula yields decent to good results for various systems 
with a minimum of effort. Even in the case of force functions without closed- 
form potential energies or discontinuous force functions approximate frequencies 
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can be easily obtained in an analytically closed form. It was further shown that 
the method can be extended to asymmetric force functions. 
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A Equivalence between the simple frequency ap¬ 
proximation formula and a version of He’s 
amplitude frequency formulation 


In He’s amplitude frequency formulation m an approximate frequency for a 
nonlinear oscillator x + f(x) = 0, x(0) = A, x(0) = 0 is given by 


Riiti) ~ Riiti) 


(26) 


where the residuals are given by Ri^(t) = (Acosfai^t)) + /(Acos^i^i). He’s 
original version used the parameters u >2 = 1, ui\ = u and t\ = £2 =0. 

In [8] Geng and Cai obtained improved results by choosing the same trial 
frequencies yet, without further justification, different times 


ii,2 = Ti j2 /12 


(27) 


where T ± t 2 = 2tt/oj1,2 which is equivalent to uii^ti ,2 = 7t/ 6. Inserting these 
values into [261 yields u\ F = ) +/(Acos(i))-^(-Aco S (f )+/(Acos( f )) 


-w 2 Acos(f )+/(Acos(f))-(-Acos(f)+/(^4cos(f)) 


a _1 

—- coinciding with J9|). This explains the hitherto only phenomenologically 


&A 


observed effectiveness of the particular choice m for the times t± t 2 - 


10 





